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Much progress has been made during recent years in understanding the 
liquid state. Most of it has been based, quite properly, upon the simplest 
possible model, that of spherically symmetrical molecules, of which the 
monatomic noble gases furnish the examples. However, most of the 
molecules with which we have to deal are polyatomic, and we are faced 
with the question of their conformity to the monatomic model. 

The structures of several monatomic liquids have been revealed by x-ray 
diffraction. This structure is expressed by a radial distribution function 
which describes the short range order around a given molecule. With it, 
one can calculate the number of molecular centers within a sphere of any 
radius, r, around the reference molecule. We! have shown experimentally 
that this is a general function of the degree of expansion of the liquid over 
its close-packed structure, and have shown, further,” how it could be used to 
express the relation between the energy of vaporization of a liquid and the 
intermolecular potential energy by means of a continuous integration, 
analogous to the summation over all lattice pairs in the case of crystals. 
The formula is, 
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where E, and E, refer to the potential energies per mole in the gas and liquid, 
respectively, V is the Avogadro number, v the molal volume of the liquid, 
W the experimental distribution function, « the potential energy between 
pairs of molecules at distances, 7. If we neglect repulsion and express e by 
—k/r*®, in accord with the theory of London,’ we can write the approximate 
equations, 
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Since the integral changes but little with temperature, we can write, 
VAEyap. ~ V*(OE/OV) ~ a, 


a constant corresponding to the van der Waals’ a, though not necessarily 
having the same numerical value as that computed from critical data. It 
is a remarkably good one, for experimental values determined from AE 
and OF /OVv agree within a few per cent. 

This approach has been extended to give a formula for the energy of 
mixing two liquids whose molecules, while possessing different radii and 
fields of force, are nevertheless mixed with maximum randomness by 
thermal agitation, a kind of solution we have called “‘regular.’’* The 
fundamental expression for the energy of mixing m, and m2 molecules in 
volume v is, (24N?/v) [ni f(11) + n3S (22) + 2mn2S (12) ], where the 
integrals are those corresponding to the first equation, above, but with the 
integrals referring to like and unlike pairs, respectively. An approximate 
solution of this equation gives, for the partial molal energy of solution of 
component 1, 


B = ¢.vi[(Api/vi)* — (Aes/ve)'"], 


where ¢ denotes volume fraction, and AE, and Ak; refer to the pure liquid 
components. For regular solutions, we have, further, the possibility of 
calculating solubility relations from £;. I wish to mention that Dr. 
Scatchard’ has contributed significantly to these developments. : 

This formula has proved far more useful than one had any right to expect 
in view of the approximations that had to be made to put it into practical 
form. It not only neglects repulsions and possible dipole interaction, but 
was derived on the assumption of spherical molecules with similar, radial 
force fields. Nevertheless, it applies remarkably well to polyatomic 
molecules, provided that they are not too unsymmetrical or different in size. 
The following illustrations are selected from many. 

Vold® measured the heat of mixing CCl, and SiCl, and found 32.8 cal. 
per mole for a 50 mole per cent mixture, while the value calculated from 
the equation was 28.0 cal. 

We predicted that I, and CCl, should give two liquid phases above the 
melting point of iodine, with a consolute temperature between 150°C. and 
170°C. with the mole per cent of iodine about 70. Experiment gave 
160.5°C. and 68 mole per cent. 

The substances SiCl, and SnI, differ so greatly as to form two liquid 
phases above the melting point of the latter. From the measured solu- 
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bility’ of SnIy in SiCl, at 59°C. we calculate 117.6 cal./cc. for Ae/v of SnIy 
while the heat of vaporization gives 118.3 cal./cc. 

But in spite of such successes, there is reason to believe that the field of 
force around large, polyatomic molecules is not best described as radial. 
The London forces are very short range, and it is the peripheral atoms 
which are most important. OsFs, for example, is very volatile although Os 
is extremely non-volatile. I have sought to throw light upon this question 
by a rather simple, direct test; a comparison of the entropy of vaporization 
of two liquids under different but corresponding conditions. My col- 
league, Dr. Pitzer,* has pointed out that two liquids which obey the theory 
of corresponding states should have equal entropies of vaporization at 
equal ratios of v,/v:, but that this presupposes, among other conditions, 
radial force fields of the same shape. 

Now there are also two other rules regarding the entropy of vaporization. 
Trouton’s rule, the historic one, states that all normal liquids have the same 
entropy of evaporation at their boiling points. But this is not strictly true, 
for there is a gradual increase with increasing boiling point, and I proposed 
years ago, as a substitute, a rule that has come to be known as the ‘“‘Hilde- 
brand Rule’’,’ according to which the comparison should be made at equal 
vapor volumes instead of pressures. 

Now two van der Waals’ liquids should have the same entropy of vapori- 
zation not necessarily at equal v,/v, but at equal v,/(v: — 0), as can be 
seen from the following equations: 

(OP/OT), = R/(v — b) = (s/dv)r; 

s, — s = R [In v, — In(v, — 9)]. 
This follows also from a.simple statistical analysis.'° We can see that for 
the noble gases there is high probability that v, is proportional to v; — 3}, 
but this is not necessarily the case for larger, polyatomic molecules, where 
the peripheral atoms are the main centers of attraction. In order to bring 
to light any significant difference between evaporation at equal v, and 
equal v,/v,, it is necessary to select liquids which differ considerably in 
liquid volume but which are otherwise as nearly alike as possible. Table 1 
shows the comparison between two pairs, first, chlorine and carbon tetra- 
chloride; second, ethane and diisopropyl. Good thermodynamic data are 
available for our purpose. 





TABLE 1 
ENTROPY OF VAPORIZATION 
As AT SAME 

REFER, alii 
LIQUID me As Vo Vo/vi T/Te FOR 
Ch 205 25.3 24.7 27.2 29.7. CCk 
Cl. 223 22.6 21 7 24.2 26.3 CCL 
C.He 148 - 25.8 25.5 28.3 31.8 CeHu* 

22.7 25.5 28.3 CoH 


C.He 160 23,2 
* Di-isopropyl, ? 
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It can be seen that the entropy of evaporation, As, for CCl, agrees best 
with that for Cl, when the comparison is made at equal vapor volumes; 
less well at equal ratios of vapor to liquid, and much more poorly at equal 
corresponding temperatures.'' The agreement with the Hildebrand rule 
is even closer with ethane and diisopropyl. 

It must be concluded, I believe, that the forces between polyatomic 
molecules are chiefly those between their neighboring parts. Although the 
assumption of similar radial fields does remarkably well for molecules of 
moderate size, it will prove increasingly inadequate as one considers larger 
and larger compact molecules. The localization to which we are accus- 
tomed in the case of coulombic forces where they exist in large molecules 
is even more necessary when considering van der Waals’ forces, in view of 
their much shorter range.'” 
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URETHANE (ETHYL CARBAMATE) THERAPY IN SPONTA- 
NEOUS LEUKEMIAS IN MICE*+ 


By L. W. Law 
RoscoE B. JACKSON MEMORIAL LABORATORY, BAR HARBOR, ME. 


Communicated May 10, 1947. 


Palliative effects due to urethane therapy have been reported in human 
cases.’ * The response to urethane in lymphoid leukemia was less pro- 
nounced and more variable than in myeloid leukemia. On the other hand 
the responses of transplantable animal leukemias (rat and mouse) were 
found by the same authors* to be very much less striking than in humans, 
if evident at all. Murphy and Sturm‘ have shown a marked effect of 
urethane on the percentage of “takes” of a transplantable lymphoid leu- 
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kemia in the rat. Recently, it has been shown in a small series of mice 
inoculated with a chronic myeloid chloroleukemia that certain palliative 
effects similar to x-ray therapy were obtained.> The effect on prolongation 
of life was not determined in any of the above experiments. 

Preliminary work is reported here relating to the effect of urethane on 
spontaneous leukemias arising in the C;s and RIL, subline B, inbred leu- 
kemic strains of mice. In this laboratory approximately 90% of mice of 
both sexes develop leukemia beginning at 6 months of age in the C;s strain. 
The majority of these leukemias are lymphoid (lymphocytic) although 
occasionally myeloid and ‘‘stem-cell” leukemias appear. Early symptoms 
of leukemia in mice of this strain can be detected by periodic palpation of 
subcutaneous lymph nodes. Usually a single lymph node is initially in- 
volved.® In all such cases diagnosed in this manner as leukemia, immature 
leukocytes have been found in the circulating blood and the animal has 
developed a generalized, systemic disease. Thus, life expectancy within 
this strain may be determined easily and a moderate degree of accuracy of 
the effect of chemotherapeutic agents on length of life of leukemic mice may 
be expected. An incidence of nearly 80% leukemia has been observed in 
the RIL strain of mice. The incidence is slightly higher among females. 
Myeloid, lymphoid and ‘‘stem-cell’” leukemias have been found and the 
majority of leukemias involve initially and principally the thymus. In 
these cases the first symptom observed is dyspnea and the animals are at 
this time in the terminal stages of the disease. 

Thirty-eight cases of spontaneous leukemia have been observed in this 
preliminary study, the majority arising in the Css strain. 

Urethane’ was administered intraperitoneally in aqueous solution. The 
daily dosage was 0.5 mg./gram of body weight. The dose per gram of 
mouse was contained in 0.0083 cc. distilled water. This sub-anesthetic 
dose did not materially effect the body weight. 

The majority of spontaneous cases reported herein received urethane 
therapy beginning from 3 to 8 days following discovery of symptoms by 
palpation. Thus, in many of the mice in the experimental series, therapy 
was instituted after the disease was well advanced. The effect of ure- 
thane was similar in all spontaneous cases studied. After 3 daily doses 
(approximately 45 mg. urethane) white blood-cell counts dropped mark- 
edly. In most cases white blood-cells counts were down to normal levels 
within 7 days and with continued therapy the counts remained within 
normal limits or the animals developed a definite leukopenia. Only 2 cases 
were observed in which during the terminal stages of the disease (day of 
death or day preceding death) the white blood-cell count returned to ex- 
tremely high levels. 

Differential whitt-cell counts obtained periodically following treatment 
of spontaneous cases showed a drop in percentage of immature cells. In 
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TABLE 1—TypicaL BLoop PICTURES IN SPONTANEOUS 
FOLLOWING URETHANE TREATMENT 
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9 Crs 38T 8 4 43,375 12 46 41 13.5) 14,800 24 40 30 11.5} 12,500 38 
9 Css 80307 [14 7 144,750 66 5 28 14.5) 36,250 57 11 30 11.5) 14,250 60 
RIL 28 8 8 29,000 32 47 21 12.0} 11,625 8 72 20 11.5] 17,375 12 
9 Css 78977 [15 4 21,750 42 24 34 11.0} 34,875 47 18 34 11.0} 14,600 44 
o’'Css 79508 [15 6 24,250 37 20 42 11.0 3,000 36 30 30 11.0 5,000 24 
9 Css 62 7 5 36,625 70 17 13 15.5 6,750 55 38 6 14.5) 14,375 58 
Ces 71 8 5 13,875 44 37 19 13.0 3,850 31 56 8 13.0 8,250 16 
Averages} 44,803 43.3 28 28.3 13.3) 15,878 37 38 22.6 12.0] 12,350 36 
9 Cos 15 11} None 14,000 42 28 30 15.5} 26,375 39 26 32 8.0} 58,875 26 
o'Css 40 | 9} None 38,250 53 23 #19 15.0} 32,375 45 9 41 13.0} 26,000 * 44 
RIL 17 10} None 68,000 6 28 64 11.0} 120,000 3 36 60 9.5 
9 Ces 79 6| None 38,725 55 1l 33 36,500 34 11 53 
9 Cos 671 6| None | 17,875 34 49 17 7,875 18 72 8 12.5) 36,500 13 








* First reading following urethane therapy at 3 days. Subsequent readings at approximately 
7-day intervals. * 


most cases this was a gradual decrease with continued therapy resulting in 
levels of only a few immature cells at approximately 9 to 10 days. These 
levels usually remained low until death. Two cases of myeloid leukemia 
in the C;s strain responded in the same manner. The amount of reduction 
was found not to be proportional to the pre-treatment immature cell 
count. 

In addition to the marked decrease in total white blood cells and im- 
mature forms there resulted also a conspicuous decrease in circulating 
lymphocytes and a corresponding increase in polymorphonuclear leuko- 
cytes. This phenomenon was described by Hawkins and Murphy’ in the 
blood of rabbits subjected to urethane anesthesia. 

Following is a typical example of response of a spontaneous lymphoid 
leukemia to the drug: Css co’ 79508 was administered urethane, 15 mg. 
daily, five days after appearance of symptoms. The white-cell count had 
fallen to 12% of the original figure after 7 days of treatment. This leuko- 
penic level was maintained by therapy throughout the life of the animal. 
The immature cells showed a reduction in absolute count to 25% of the 
original after 13 days and to 4.8% (2% immature cells) of the original after 
18 days of treatment. This level was maintained by therapy. The typi- 
cal decrease in circulating lymphocytes and increase in polymorphonuclear 
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LEUKEMIC MICE OF THE Css AND RIL STRAINS 
(SuB-ANESTHETIC Dose) 





FOLLOWING URETHANE THERAPY* 
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36 15 9.5} 8,000 21 56 23 12.5) 9,625 22 53 22 10.5/8625 41 38 21 8.0 

65 11 9.3} 7,500 31 67 2 126:8)4950..8 & 2 8.0 

35 «6 14.5] 11,825 38 58 4 11.5/80,750 25 71 4 77 

66 12 11.5] 6,000 16 82 2 SED 

49 11.8 10.5] 8,046 27 63.3 9.6 9.9/14,095 19.466 12.6 8.4 
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33 54 10.5 














t Myeloid leukemia. All others lymphoid. 
} Urethane therapy started in terminal stages 12 days after appearance of symptoms. 


leukocytes was evident after 7 days of treatment. After 27 days the circu- 
lating lymphocytes had fallen from 37 to 8% and the polymorphonuclear 
leukocytes had increased from 20 to 84%. 

In contrast to the results obtained in human cases! there was not ob- 
served a general rise in hemoglobin levels. In 2 treated mice only was 
there observed a rise in hemoglobin levels. In 2 Css 78977 the hemo- 
globin value rose by an average of 10.7% (100% = 14.5 g. Hb per 100 cc.) 
after 19 days of therapy. In o Css 81 the hemoglobin value rose by an 
average of 11.1% after 17 days of treatment. Indications are, however, 
that the fall in hemoglobin levels in urethane treated leukemias is not so 
precipitous as in the untreated controls. This is more evident in spon- 
taneous cases given urethane immediately upon discovery of symptoms. 
These cases are not included in this series in which hemoglobin levels were 
relatively low when therapy was initiated. (See table 1 for typical re- 
sponse of blood counts in spontaneous cases of leukemia.) 

Marked changes in the blood picture were evident 24 hours after ad- 
ministration of urethane (sub-anesthetic dose). The leukemic blood pic- 
ture returned immediately upon cessation of treatment. 

In mice given urethane within 3 to 5 days after discovery of palpable 
nodes, the subcutaneous node (or nodes) regressed completely within 2 to 
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3 days of continued therapy. In cases where urethane was given 8 days or 
later after appearance of symptoms a definite decrease in size resulted. 
The difference in mean weights of subcutaneous lymph nodes (axillary, 
inguinal and cervical) in the control and urethane-treated cases was 248.8 
mg. where ¢ = 2.7 and P = 0.02%. 

There was a definite decrease in the size of the spleen in spontaneous cases 
already showing splenic involvement upon therapy. In cases where treat- 
ment was begun before the spleen was palpable the spleen did not become 
infiltrated greatly. The difference in mean weights of the spleen in control 
and urethane mice obtained at autopsy was 324.0 mg. where ¢ = 1.29 and 
P < 0.3%. There was considerable variation in splenic weights in the 
urethane-treated series probably as a result of institution of therapy at 
various time intervals after appearance of symptoms. 

TABLE 2 


COMPARISON OF MEAN WEIGHTS OF SPLEEN, LYMPH NODES AND THYMUS IN URETHANE- 
TREATED AND CONTROL SPONTANEOUS LEUKEMIAS 


- —————-WEIGHTS IN MILLIGRAMS————————. 
NO. SPLEEN LYMPH NODES* THYMUST 
Uretharie-treated ‘i 731.4 194.3 73.7 
Controls 13 1055.4 479.1 530.8 
Difference (means) = 324.0 284.8 457.1 
t= 1.29 2.7 3.7 
Pp = < 0.3% 0.02% . < 0.01% 


* Including all axillary, inguinal and cervical lymph nodes. All weights obtained at 
death of animals. 

+ Thymus weights on 8 experimental and 13 control animals were obtained principally 
from RIL, subline B, pedigreed mice and are unrelated to spleen and lymph node 
weights obtained principally from C;s pedigreed mice with systemic disease originating 
in subcutaneous lymph nodes. 


The mean length of life of untreated spontaneous leukemias in the Css 
strain was 21.8 + 7.15 days and in urethane-treated leukemias 35.9 = 8.76 
days following discovery of symptoms by palpation. In one case, 2 Css 
80302 life was maintained for 52 days. Although the difference of the 
means is not significant in the series reported here, indications are that if 
urethane is given early, following discovery of symptoms, a greater life ex- 
pectancy may be obtained. 

Leukemic mice from both the C;s and RIL strains (8 animals), in which 
there was initial thymic involvement resulting in dyspnea followed usually 
by a rapid, generalized systemic course of the disease, were given urethane 
treatment immediately upon discovery of the dyspnea. There resulted a 
marked depression in the white-cell count, a decrease of immatures in the 
circulating blood and in most cases a sudden disappearance of dyspnea. 
In these mice there was a significant increase in the length of life, 22.7 + 
3.1 days compared with 7.7 + 5.34 days, for untreated controls. At 
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autopsy it was found that the thymus had regressed, was soft, spongy and 
yellowish-brown in color. The difference in mean thymus weights between 
the urethane-treated and control series was 457.1 mg. where ¢ = 3.7 and 
P < 0.01% (see table 2). 

The following histological picture obtains in leukemic mice receiving 
daily sub-anesthetic doses of urethane: In the lymph nodes of mice receiv- 
ing urethane shortly after appearance of symptoms there was observed in 
most cases an intense necrobiosis with pyknosis, karyorrhexis and chroma- 
tolysis. This was usually generalized throughout the nodes. However, 
the process was localized in leukemic mice in which therapy was begun in 
the terminal stages of the disease. Where therapy had been continued for 
a relatively long period of time there were found areas of necrosis in the 
lymph nodes. In several cases moderate to intense reticulosis and fibrosis 
were observed. The reactions of splenic tissue to urethane were similar but 
not so intense as in lymph nodes. 

The most pronounced effect of the drug has been observed in the thymus 
where intense generalized necrobiosis, large areas of necrosis and fibrosis 
and severe hyperemia were observed. 

Phagocytosis was very prominent in the livers of treated mice. Large 
areas of perivascular fibrosis in regions of leukemic infiltration were pres- 
ent. In some hepatic cells there was a moderate degree of vacuolar de- 
generation. Hyperemia was pronounced in some livers. 

In o& RIL 28 in which there was severe thymic involvement and in which 
urethane therapy was started 8 days after appearance of symptoms there 
occurred moderate leukemic infiltration in the lung. There was observed 
in these areas of infiltration a very intense necrobiosis. 

More complete details concerning the histology of urethane-treated 
tissues will be given later. 

Summary.—Definite palliative effects of daily sub-anesthetic doses of 
urethane have been observed in 18 cases of spontaneous leukemia in the 
Css and RIL inbred strains of mice. The majority of these generalized 
systemic leukemias were lymphoid although the response of two myeloid 
leukemias was similar. The effects produced were: (1) a pronounced fall 
in total white-cell count to or below normal limits and a maintenance of 
these levels with continued urethane therapy; (2) a marked reduction in 
the number of immature cells in the circulating blood; (3) a possible sta- 
bilizing effect on the hemoglobin levels; and (4) a pronounced diminution 
in the size of enlarged subcutaneous lymph nodes and spleen. The life 
expectancy of urethane-treated leukemics is greater than control leukemics 
but not of statistical significance. 

In a series of 8 leukemic mice of the RIL strain in which there was initial 
and principal thymic involvement the same palliative effects were ob- 
tained. In addition there was relief from symptoms of dyspnea and tho- 
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racic enlargement, a statistically significant decrease in the size of the 
thymus, and a statistically significant lengthening of life. 

A generalized necrobiosis in the spleen, thymus and lymph nodes followed 
by localized necrosis in animals receiving long continued therapy has been 
observed. Moderate to intense reticulosis and fibrosis were evident in 
cases receiving treatment early and over a relatively long period of time. 


* This work has been aided by grants to the Roscoe B. Jackson Memorial Laboratory 
by the Jane Coffin Childs Memorial Fund and the National Advisory Cancer Council. 

Tt Technical assistance of Mr. Lester E. Bunker, Jr., and Miss Betty-Ann Norris is 
gratefully acknowledged. 
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INFLUENCE OF HUMIDITY ON THE SURVIVAL OF DIFFERENT 
CHROMOSOMAL IT YPES IN DROSOPHILA PSEUDOOBSCURA 


By M. J. Heuts* 


UNIVERSITY OF LOUVAIN, BELGIUM, AND DEPARTMENT OF ZOOLOGY, COLUMBIA UNI- 
VERSITY, NEW YORK 


Communicated May 7, 1947 


Populations of Drosophila pseudoobscura show a high variability with re- 
spect to the gene arrangement in the third chromosome. Several gene 
arrangements, which must have arisen from each other by inversions of 
blocks of genes, occur frequently in the same territory. Populations of 
different localities very often differ in the relative frequencies of the gene 
arrangements,' and in at least some localities, the frequencies change also 
from month to month, the changes being connected with tht succession of 
year’s seasons. In particular, seasonal changes are observed in the popu- 
lations of Pifion Flats and Andreas Canyon, on Mount San Jacinto in 
California. The Standard gene arrangement increases in frequency during 
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the hot period of the summer, remains rather constant during autumn and 
winter, and decreases in frequency during spring. The Chiricahua arrange- 
ment undergoes a cycle opposite in sign to that of Standard, and the Arrow- 
head arrangement shows relatively little change. These data suggest that 
carriers of different gene arrangements differ in adaptive value, and that 
the seasonal changes represent adaptive responses of the populations to 
changing conditions of the milieu. Artificial populations kept in spe- 
cially constructed ‘‘population cages’’ bear out the hypothesis of natural 
selection, in so far as, at relatively high temperatures (above 21°), changes 
are induced which seem to parallel those observed in natural habitats in 
summer. Thus, if the initial population of a population cage contains less 
than 50 per cent of Standard and more than 50 per cent of Chiricahua 
chromosomes, the former increase, and the latter decrease in frequency, 
until an equilibrium is established at about 70 per cent Standard and 30 
per cent Chiricahua. Similarly, Standard is, at high temperatures, rela- 
tively superior to Arrowhead, and Arrowhead is superior to Chiricahua. 
At temperatures below 21°, the adaptive values of the carriers of all gene 
arrangements appear to be alike.* 

The experiments just referred to describe, however, the net differences 
between adaptive values of the different: gene arrangements. It is not 
known at just what stages of the developmental cycle the differential sur- 
vival takes place. Furthermore, in the experiments so far published, the 
carriers of the Chiricahua gene arrangement have never proved superior in 
adaptive value to those with Standard and Arrowhead. Yet, in the nat- 
ural populations, Chiricahua chromosomes do increase in frequency at the 
expense of the others during the spring season, It follows that, an eco- 
logical niche should be found in which the carriers of Chiricahua will be 
relatively superior. One of the environmental agents never tried before 
which may have an effect on the survival of the chromosomal types, is 
humidity. The pupal stage is especially likely to be sensitive to humidity 
variations. Among the known facts which point in this direction, suffice it 
to mention that, differential mortality of wild type pupae and pupae 
homozygous for certain mutants has been found in Drosophila melanogaster 
at low humidities.* 

In the experiments to be described below, Drosophila pseudoobscura 
pupae, homozygous for Standard, Arrowhead, and Chiricahua gene arrange- 
ments have been used. The ancestors of all the experimental flies came 
from the Pifion Flats locality in California. At least ten strains with each 
gene arrangement were intercrossed, so that the experimental flies, though 
structurally homozygous, were genically heterozygous (the importance of 
which is discussed by Wright and Dobzhansky*). The parental flies were 
placed in population cages, some hundreds of males and females of each 
type in a different one. The cages were kept in a constant temperature 
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Ordinates: percentages of hatched pupae. Abscissae: relative humidity in per cent. 
Solid circles—Standard pupae; open circles—Chiricahua pupae; triangles—Arrowhead 
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56% 
56% 

0% 

0% 

0% 


room at 19°, and at relative humidities ranging from 45 to 50 per cent. 
When pupae began to form in the cages, they were extracted individually 
by means of a needle, and transferred into glass vials, one hundred pupae 
The vials were closed by cheese cloth held by a rubber band. 


per vial. 
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FIGURE 1 


SECOND EXPERIMENT 
FLIES 


TABLE 1 

FIRST EXPERIMENT 
PUPAE FLIES % PUPAE 
300 248 82. 400 
300 276 92. 400 
300 220 73. 400 
300 278 92.6 300 
300 233 Ce he 300 
300 256 82. 300 
500 «6419S 8838.. 300 
500 376 = = 75. 300 
300 265 88.¢ 300 
300 
si ne <p 200 
300 122 40. 700 
300 94 31.3 700 
300 215 71.0 700 
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468 
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Only very young and light pupae were taken. The vials were placed in 
desiccators at the desired humidities. Five relative humidites were used, 
namely 100%, 92%, 76%, 56% and 0%. They were obtained by placing 
on the bottom of the desiccators, distilled water, K2SO,, NaCl, NaBr in 
oversaturated solutions, and anhydrous CaCh, respectively (according to 
Ludwig and Landsman‘). The desiccators with pupae were placed at a 
constant temperature of 25°C. 

The experiments were replicated twice. The results obtained are re- 
ported in table 1, and represented graphically in figure 1. It can be seen 
clearly that, at the temperature of 25°C., Chiricahua pupae are more viable 
than Standard, and the latter more viable than Arrowhead at 100% hu- 
midity. At 92% humidity, Standard is superior to the other two. At 
76%, Standard and Arrowhead are alike, but both of them are superior to 
Chiricahua. At lower humidities, Arrowhead is superior to both Standard 
and Chiricahua. Thus, Chiricahua is most favorable, and Arrowhead 
least favorable, at 100% humidity, but the relations are reversed at 0%. 
The two series of experiments showeci the same hierarchy of hatchabilities 
of the pupae with the three gene arrangements, although in the first experi- 
ment the hatchabilities were lower in almost all humidities than they were 
in the second experiment. The difference must have been due to more 
favorable environmental conditions in the second than in the first series of 
population cages in which the pupae were reared. 

The results of the experiments just reported suggest that, flies with 
Chiricahua chromosomes are relatively better adapted to climates with 
high humidities, and Arrowhead flies to low humidities, Standard being 
intermediate. The temporary superiority of Chiricahua observed during 
the spring season at Pifion Flats might, then, be due to the relatively high 
humidities prevailing at that season. The validity of this hypothesis must, 
of course, be tested by studying the survival rates of pupae at different tem- 
peratures? 

I wish to express my gratitude to Prof. Th. Dobzhansky for his generous 
hospitality throughout the course of this investigation, and for his con- 
tinuous advice and encouragement. 


* Fellow of the Belgian American Educational Foundation, Inc. 

1 Dobzhansky, Th., and Epling, C., Carnegie Inst. of Washington Publ. 554 (1944). 
2 Dobzhansky, Th., Genetics, 28, 162-186 (19438). 

3 Wright, S., and Dobzhansky, Th., Genetics, 31, 125-156 (1946). 

4 Geisler, F. S., Amer. Natur., 76, 223-224 (1942). 

5 Ludwig, D., und Landsman, H. M., Physiol. Zool., 10, 171-243 (1937). 
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PERTURBATIONS OF DISCONTINUOUS SOLUTIONS OF NON- 
LINEAR SYSTEMS OF DIFFERENTIAL EQUATIONS 


By NorMAN LEVINSON: 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


Communicated May 3, 1947 
Theorems about solutions of the system 


e at. Sek A. ee eee | (1) 
where the X;, are regular functions of ¢ for small |e|, are classical. More 
recently studies of cases of (1) subject to the various conditions that allow 
one or several X; to become large as « —> 0 have been studied.! 

In all the studies just referred to it is assumed that the system (1) has a 
solution with a continuous derivative in case e = 0; sufficient conditions 
are given for this also to be the case fore # 0. There are cases of con- 
siderable interest where, when e = 0, the system (1) has only a discontinu- 
ous “‘solution.”’ It is known as a practical matter that in these cases when 
¢ ~ 0 but small the system may have a continuous solution which ap- 
proaches the discontinuous one as e —> 0. Practical applications of this 
have been made by the Russian school of non-linear mechanics (Mandle- 
stam, Chaikin, Lochakow). A rigorous treatment of this case appears so 
far to have been given’ only for a special relaxation oscillation problem in- 
volving two unknown functions (so that the simple topology of the phase 
plane could be exploited). 

Here the following system of arbitrary order is considered: 


dx; du ‘ 

a IG toot = 12m 

du du 

oP + ou +h= 0. (2) 


In (2), fi, os, g and / are all functions of x, ...x,, u, and ¢ which we shall 
denote hereafter by f;(x, u, /), etc. We shall first show that (2) is general 
enough to include certain cases of obvious interest: 
Example 1. The van der Pol equation, with a change of time scale, be- 
comes 
d*u 


€ 


dt? 





+ — + = 0; 
U it u = VU; 


Example 2. The Rayleigh type equation, 
ea + (xe? -1l)e+x4+x'3 = 0, 
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is not included in (2). However if we differentiate it and set = u we get 


dx 
— = %, 


dt 
oH + ut — 1) 4 uw + Bx%u = 0 
it? ou dt Uu xu = 


which is of the form (2); 
Example 3. Consider the system* 


d. 

2 = H,(x, w, t),4 = 1, ..., ®, 
d 

c= G(x, w, t). 


In case the right members above are not linear in w, this system can be 
brought to the form (2) simply by differentiating the last equation only 
with respect to ¢. 
In case the right members are linear in w, we can introduce u defined by 
= wand the system assumes the form 


ise, 
Ot fle, oe + $:(x, 0), 
2. 
— = ee, “+ h(x, t), 
‘ “it? 


which is a special case of (2). 
To formulate our result precisely we consider along with (2) where « > 0, 
the degenerate system 


dy; 
Ot = fly, 04) S | + by, 8,4 = 1,2, 


d 
By, », D7, + ho, 0, f) = 0, (3) 


which is (2) with e = 0. We shall regard a solution of (3) as a curve Cp in 
E,,+2, the n + 2 dimensional space (y, v, ¢). We shall assume that there is 
an open continuum D in E,,+2 which contains the solution Cy and in which 
fi %y g and h are continuous and have continuous first order partial de- 
rivatives with respect to y;, v and ¢. 

Definition. Then + 1 functions y,(¢), v(t), denoted briefly by the n + 1 
dimensional vector Y(t), is said to be a solution of (2) the degenerate system 
(3) in the interval a S$ ¢.S 6 if: 
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(1.1) Y(é) is continuous and possesses a continuous derivative except 
possibly at a finite number of interior points, 7,, of the closed interval (a, 8) 
and except at the 7;, Y(t) = (y,(t), v(t)) satisfies (3). 

(1.2) Except at the points 7;, 

g(y(t), v(t), t) # 0. 

(1.3) Both Y(7; — 0) and Y(r; + 0) exist, and we denote them by Y- 
and Y+, respectively. Similarly we define y-, v-, etc. By f;- we mean 
id=; v-, Tj), etc. 

(1.4) For each j, g— = 0; 

(1.5) For each j, g+ # 0; 

(1.6) Denoting by y(v, 7;) the solution, y;, of 

dy; 


= i( 1, 0, Ts 
dv f 9 Tj) 


for which y(v-, 7;) = y- and y(v+, 7;) = y+, let 
S g(y(v, 7;), v, 7;)dv = 0. 
(1.7) There exists no v interior to the interval (v-, v+) such that the 


integral in (1.6) vanishes if v+ is replaced by % in the upper limit of the 
integral. . 


Og”, =<, Of ; 
(1.8) Let J denote a +> > fi. Replacing y by y-, v by v— and t by 7;, 
v i=1 OY; 
I becomes J—~. We assume that 
(7) > ©: 


(1.9) For small 6 > 0, 
g(x(7; — 6), v(7; = 5), | ies 6) > 0. 


The curve Cy in E,,+2 is determined by Y(#) for ¢ ¥ 7; and by y(v, 7;) 
where v goes from v— to v+ when ¢ = 7;. Thus Cp is continuous, and Cp 
possesses a tangent except at the several points Y— and Y+. 

The above rather complicated definition of a discontinuous solution of 
(1.3) is justified by the following result. 

THEOREM I. Let the degenerate system (3) have a solution Co, fora StS 
B. Let X(t) = (x(t), u(t)) bea solution of (2) for « > 0 such that 





X(a) — Y(a)| = Llx(a) — yi(a)| + \u(a) — v(a)| < 6 


du(a) dv(a)| Ss bo 
dt dt | 


€ 
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If €, 6; and 5, are sufficiently small, X(t) is a solution of (2) overa St S 8B. 
Moreover as €, 6; and 5, — 0, the curve representing X (t) in E,,42 approaches Cp. 


du(8) _, dv(8) 


In particular X(8) > Y(B). Also : 
dt dt 


Thus we see that the existence of a solution, even though discontinuous, 
of the degenerate system (3) implies the existence of nearby solutions of the 
system (2). 

Let us denote the initial values of one of the » + 1 codrdinates of X(a) 
and Y(a) by a. Then we can regard the values of X and Y att = Bas 
functions of Banda. We have 

THEOREM II. Ase, 6; and & — 0, 


dX(B, a) _, 9¥(B, a) 
oa a 





Also as €, 51, 55 — 0, 
O du(B, a) _ 0 do(B, a) 
Oa dt Od «dt 


: du 
Moreover denoting the value of i“. at a by c, we have as e— 0 


OX(Be) 4, 2 dulBe) _, 5 
Oc Oc dt 

In case the degenerate system (3) has a periodic solution and the Jacobian 
associated with the determination of this solution by varying initial condi- 
tions is distinct from zero, then it follows from Theorem II that the corre- 
sponding Jacobian for (2) also will not vanish and therefore (2) will also 
have a periodic solution.‘ 

The treatment can be generalized readily to the case where not one but 
several or even all the equations of (2) assume the same form as the last 
(nm + 1st) equation and moreover where the terms f, ¢, g and h depend on e€ 
(but approach finite limiting values as « — 0). 

A further generalization is the following system: 


dx; 


= 3 fale, wt ae y+ be ub €), 4 eat RR 
dt j=l 


& 
at Leute ut, oe | + hile, ut, 6) = 0, t= 1,...m. 


Let us denote the square matrix 


(gis(x, u, t, €)) 








218 PHYSICS: G. POLYA Proc. N. A. S. 


by G and its determinant by |G). We assume that the degenerate system 
(« = 0) has a solution which can be continued up to ¢ = 7 where 7 is such 
that for the solution in question |G(x(r), u(r), 7, 0)| = 0. Part of our 
sufficient condition for the perturbed system to have a solution is the re- 
quirement that the characteristic equation 


\gas(x(r), u(r), Ty 0) ce bij] sad 0 
have \ = 0 asa simple root. 


1M. Nagumo, Uber das Verhalten der Integrals von Ay” + f(x, y, y’, ) = 0 fiir 
X40. Proc. Phys. Math. Soc. Japan, 21, 529-534 (1939). I. M. Volk, A Generaliza- 
tion of the Method of Small Parameter in the Theory of Non-Linear Oscillations of Non- 
Autonomous Systems. C. R. (Doklady) Acad. Sci. U.S.S.R., 51, 437-440 (1946). Volk 
considers (1) where the X; are meromorphic functions of ¢ for small ¢ and periodic in ¢. 
K. O. Friedrichs and W. R. Wasow, Singular Perturbations of Non-Linear Oscillations. 
Duke Math. Jour., 13, 367-381 (1946). Here the X; are not functions of ¢ and for 7 S 
n — 1are not functions of «. X, contains ¢ in the form of a factor 1/e. 

2D. A. Flanders and J. J. Stoker, The Limit Case of Relaxation Oscillations, Studies in 
the Linear Vibration Theory, New York Univ., 1946. 

3 This is the system, except that ¢ is not necessarily excluded from the right members, 
which is considered by Friedrichs and Wasow, loc. cit. 

4 See Friedrichs and Wasow, and Volk, Joc. cit, for continuous cases where right 
members do not and do, respectively, depend on ¢. 


A MINIMUM PROBLEM ABOUT THE MOTION OF A SOLID 
THROUGH A FLUID 


By G. POLya 
STANFORD UNIVERSITY 
Communicated April 25, 1947 


1. An incompressible frictionless fluid of uniform density p fills the 
whole space outside a moving solid and is at rest at infinite distance. 
The motion of the solid is one of pure translation. The magnitude of the 
velocity is U, its direction cosines with respect to a codérdinate system fixed 
in the solid \, u,v. The kinetic energy of the fluid is of the form 


T = 1/2MU’. 


The quantity M, called the virtual mass, depends on the direction of the 
velocity : 


M/p = Ad? + Bu? + Cv? + 2A’pv + 2B’rr + 2C’An. 


A, B, C, A’, B’, C’ are uniquely determined if the shape and size of the solid 
and the relative location of the coérdinate system and the solid are given. 
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A closer study of the dependence of A, B, C, A’, B’ and C’ on geometric 
data may seem desirable.' Taking a first step in such a study, we consider 
the average virtual mass M, obtained by averaging M over all directions , 
u,v and assuming p = 1: ~ 


M =(A+B+0)/3. 


M is independent of the location of the codrdinate system and depends only 

_on the size and shape of the solid. It is easy to show that of all ellipsoids 
with given volume the sphere has the minimum average virtual mass. It would 
be natural to suspect that this statement remains true if for ‘‘ellipsoids”’ 
we substitute “‘solids.’’ At any rate, I shall prove the analogous general 
theorem in two dimensions. 

2. We consider now the two-dimensional motion of an incompressible 
frictionless fluid of uniform density p that fills the space around a cylinder 
of infinite length. The motion is parallel to a plane, the plane of the com- 
plex variable z, that is perpendicular to the cylinder and intersects it in a 
closed curve C (the notation of section 1 has been dropped). The exterior 
of Cis mapped conformally onto the exterior of the unit circle in the ¢-plane 
so that the points at infinity correspond. Thus, z moving outside C is 
represented by the series 


es Ses ae ) (1) 


convergent for |¢/ > 1. The number ) is positive. 

We begin with the case in which the motion of the fluid at infinite 
distance is parallel to the real axis and has the velocity U (uniform flow 
disturbed by a fixed cylindrical obstacle). The corresponding motion in 
the ¢-plane, around a circular cylinder and with velocity Ud at infinity, 
has the complex potential 


1 
x’ = urs + :). (2) 


Yet (2) represents also the complex potential for the z-plane provided that 
z and ¢ are linked by the mapping (1) that transforms streamlines into 
streamlines and, especially, the unit circle of the ¢-plane into C. 

3. To the motion just considered we add a uniform velocity U directed 
along the negative real axis. We obtain thus a new motion (disturbance of 
a fluid which is at rest at infinite distance by a cylinder moving through it 
sidewise to the left). The complex potential of this motion is obviously 


Cae Me eae OS u|r(s os 7). | (3) 


where z and ¢ remain linked by (1). (Of course the cofrdinate system re- 
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mains fixed with respect to the solid.) The velocity at the point z is @, 
conjugate to 
—~' be (4) 
i | 


The kinetic energy of a layer of the fluid, of unit thickness and parallel to 
the z = x + ty plane, is 


1 roy 1 
Lo ff tut dx dy = 5 MU*?. (5) 


The integral is extended over the exterior of C, and M is the virtual mass 
per unit height. From (4) and (5) we obtain 


: 1 \d¢ r 
M/p = A 1 — — — 1 dxdy (6) 
CJ dz 


dz LV 
re - a(1 — ) dé dn; 
dg i 


the latter integral is extended over the exterior of the unit circle in the ¢ = 
+ in plane. Introducing (1) and polar coérdinates, we obtain from (6) 
in the usual way that 


d 


w= 
dz 


~ 

ll 

oo 
= 

al 
oa 

| 
~~ 
eee” 
aye 
ti) ~~ 


II 


M/p = rX(\o. — 1/2 + Ne! + 3lcs|2 + ...). (7) 


Now the area of C or, what is numerically the same, the volume V of the 
moving cylinder per unit height is 


V = wd*(1 — Je]? — 2leo|? — 3losi2 — ...). (8) 


' 


This is well known and obtained by a computation analogous to the one 
just sketched. It follows from (7) and (8) that 


V + M/p = 2rd7(1 — Rey). (9) 


where &c, denotes the real part of c. 

4. Now, we wish to obtain M,, the virtual mass per unit height 
corresponding to a direction of the velocity that includes the angle a with 
the direction just considered. We reduce this problem to the foregoing by 
a rotation, introducing the new complex variables 2’ and ¢’, 


, ta se! 


2 =e 2, t’ = e%. 


We obtain from (1) that 





2ia sia 

ce C2€ 

2! = M(¥ + coe** + ee + a +. ae rw) ) (1’) 
‘ths 
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Substituting c,e”** for c, in (9), we obtain 
V + M./p = 2ed2(1 — Ree**) (9’) 
and hence 
V+ Mate/2/p = 2ed2(1 + Rere"**) (10) 
We define M, the average virtual mass per unit height by” 
M = ('/emp) fo" Mada = ('/2p)(Ma + Ma+n/2)- (11) 


(M has, in fact, the dimension of an area, and so has V.) From (9’), (10) 
and (11) we find finally 


V + M = 2d? (12) 


5. Now dA is the so-called outer radius of C (that is the radius of the 
circle onto the exterior of which the exterior of C is so mapped that the 
points at infinity correspond to each other with unit magnification). It 
follows from (8) (and is well known) that 


V < mr’. 
unless Cisacircle. Therefore, by (12), 
M>vV 


with the same proviso. For the circle, however, M = V. Thus, we have 
proved that of all cylinders having the same area of the cross-section, the 
circular cylinder has the minimum average virtual mass per unit height. 

We can derive another result from (12): the average virtual mass per 
unit height decreases by symmetrization. Indeed, we know that symmetriza- 
tion leaves V unchanged and decreases the outer radius X.” 


1 This is suggested by a systematic study of the dependence of the capacity on geo- 
metric data which has been undertaken recently by Mr. G. Szegé and the author. 

2 See G. Pélya and G. Szegé, ‘‘Inequalities for the Capacity of a Condenser,” Amer. 
Jour. Math., 67, 1-32 (1945), especially pp. 13-14. 
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